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We investigate the canonical quantization of the electromagnetic field on the Kerr 
background. We give new expressions for the expectation value of the electromag- 
| netic stress-energy tensor in various vacua states and give a physical interpretation 

■ of the separate terms appearing in them. We numerically calculate the luminosity in 

these states. We also study the form of the renormalized stress-energy tensor close 
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to the horizon when the electromagnetic field is in the past Boulware state. 



INTRODUCTION 



o . 

VO . The properties of quantum states of physical interest in the Kerr background has been 



broadly investigated by studying the scalar field. The study of the quantized electromag- 
netic field in this background has been studied much less (notably, in [1], [2] and [3]). In 
this paper we perform the canonical quantization of the electromagnetic field in the Kerr 



background and calculate the renormalized expectation value of the stress-energy tensor 
(abbreviated as RSET) when the electromagnetic field is in various states of physical in- 
terest. We analytically calculate new expressions which are symmetric under parity for the 
expectation value of the stress-energy tensor in these physical states. These expressions cor- 
rect those given in Candelas, Chrzanowski and Howard [1], hereafter referred to as CCH, 
which are not symmetric under the parity operation even when the state is and even though 
the parity operation is a symmetry of both the Kerr metric and Maxwell's equations. We 
also numerically calculate differences between two states of the RSET. We show that the 
difference in the RSET between a Hartle-Hawking-type state and the past Boulware state 
is exactly (minus) thermal near the horizon, unlike a result in CCH. We also study the rate 
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of rotation of this difference, for which there is no unanimous consensus in the literature, 
and find that it is rigidly rotating with the horizon to at least second order in the distance 
from the horizon. 

This paper is organized as follows. We outline the classical theory of the electromagnetic 
field on the Kerr background in Section II. In Section III we lay out a canonical quantization 
of the electromagnetic potential and field on this background. In the following section we give 
a description of the main physical quantum states. Section V is split into two subsections. 
In the first one we derive expressions for the expectation value of the electromagnetic stress- 
energy tensor when the field is in these states. We also show that similar expressions given 
by CCH result in RSETs which are not symmetric under the parity operation. In the second 
subsection we give a physical interpretation of the various sets of terms appearing in these 
expectation values. In Section VI we calculate the luminosity of the Kerr black hole in the 
past Boulware and past Unruh states for the spin-1 case. In the last section we study the 
behaviour close to the horizon of the difference in the RSET between a Hartle-Hawking-type 
state defined in CCH and the past Boulware state. 

The method we followed to solve the radial Teukolsky equation is described in Appendix 
A, whereas we solved the angular equation as described in [4]. Further details can be found 
in [5]. In Appendix B we extend to the spin-1 case an asymptotic analysis in [6] for the 
radial solutions close to the horizon in the scalar case. 

We use geometrized Planck units: c = G — H — 1, and follow the sign conventions of 
Misner, Thorne and Wheeler [7]. All figures have been obtained for the values: a = 0.95M 
and M — 1; the Boyer-Lindquist radius of the event horizon for such values of the black 
hole parameters is r + ~ 1.3122. 

II. CLASSICAL THEORY 

The Kerr metric corresponding to a black hole of mass M and intrinsic angular momen- 
tum a per unit mass as viewed from infinity is, in the Boyer-Lindquist co-ordinate system 
{t,r,6, </>}, given by 

ds 2 = ~ (dt - a sin 2 6&<P) 2 + ^ [(r 2 + a 2 )d0 - adt] 2 + ^dr 2 + Zd6 2 (2.1) 

ZA 
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where £ = r 2 +a 2 cos 2 9 and A = r 2 — 2Mr+a 2 . The Kerr metric has co-ordinate singularities 



at A = 0, which has for solutions r± — M ± v 'M 2 — a 2 . The hypersurface r = r + is the 
outer event horizon. The surface gravity on the outer [inner] horizon is k+[ac_], given by 
K± = {r±-r^)/{2{rl + a 2 )). 

The Kerr space-time is stationary and axially symmetric; the two Killing vectors asso- 
ciated with these two symmetries are k j = d/dt and k^ = d/dcf) respectively. We can 
construct another Killing vector as k^ = k^ + Q + k(^ with Q + = a/(ri + a 2 ). The surface 
where the Killing vector k j becomes a null vector is called the stationary limit surface. The 
vector kj is timelike outside this surface and is spacelike in the region between the event 
horizon and the stationary limit surface. This region is called the ergosphere. The surface 
where the Killing vector k f %> becomes null is called the speed-of-light surface. The vector 
kjj? is timelike between the event horizon and the speed-of-light surface, and it is spacelike 
outside it. 

An observer who moves along a world line of constant r and 9 with an angular velocity 
Q = d(f)/dt relative to the asymptotic rest frame has a tetrad {e( f ), e( r ), e^, e^)} associated 
with him. Such an observer sees no local change in the geometry and is therefore considered 
a stationary observer relative to the local geometry. If his angular velocity is zero, and 
therefore he moves along a world line of constant r, and 9, he is a static observer (SO) 
relative to radial infinity. A SO moves along the integral curves of kj. If we require and 
e(0) to be parallel to d/dr and d/d9 respectively, we then find that the two other vectors in 
the tetrad of a stationary observer are given by 



= 1 (d_ Q d_ 

e(t) v /\ gtt + 2n gt(P + n 2 g H \ \dt + 

1 1 



y/\g tt + 21^ + Wg^\ ^ - g tt g H L 



■{914 + &9h)^1 + (9tt + n 9uf>)-^7 



(2.2) 



The quantity = pk^ = p a k^ where p is the 4-momentum of a certain observer, is the 
component of angular momentum of that observer along the black hole's spin axis. This 
quantity is conserved for geodesic observers. The only stationary observers for whom this 
quantity is zero are those with an angular velocity H Z amo — ~9t<f>l ' 9<f>cj>- These observers are 
called zero angular momentum observers (ZAMO). These observers are the closest analogue 
to the static observers in Schwarzschild space-time on the Kerr-Newman space-time in the 
sense that their 4-velocity is orthogonal to k$, the hypersurfaces of constant t. The angular 
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velocity ^zamo of the ZAMOs as they approach the horizon tends to Q + , which can therefore 
be interpreted as the angular velocity of the horizon. The stationary observers whose angular 
velocity is constant and equal to fl + are observers that follow integral curves of fejf and are 
called rigidly rotating observers (RRO). The last orthonormal tetrad that is of interest to 
us is the Carter orthonormal tetrad [8], which corresponds to a stationary observer (2.2) 
with angular velocity ^carter = a/(r 2 + a 2 ). 

The Newman-Penrose (NP) formalism ( [9]) is based on four null basis vectors: 
{/, n, m, m*}. In this paper we will use the Kinnersley tetrad: 
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(2.3a) 
(2.3b) 
(2.3c) 



where p = —\j(r — ia cos 9) is one of the spin coefficients of the Kerr metric in the Kinnersley 
tetrad. The other spin coefficients and NP scalars can be found in [10]. Of interest to us 
are the NP Maxwell scalars: 0_i = Fi m , O = {Fin + Fm*m) /2 and +1 = F m * n , where 
F a p = A/3-u — A a ./3 is the anti-symmetric Maxwell tensor and A a is the electromagnetic 
potential. From these equations, one can obtain the operator Kt which maps the potential 



onto the NP scalars, i.e., 



The explicit form of the operator K% can be found in [11]. 

It is easy to see that the effect of the parity operator V : 
NP Kinnersley tetrad is 



(2.4) 



f 7r — 6, d) + it) on the 



VI = I, 



Vn = n, 



Vm = —m* 



(2.5) 



In the limit r — * +oo, the Kinnersley tetrad becomes 

1 ,„ 



I 



-e t + e r 



n 



y e t + e r ) , 



m 



V2 



.e + ie 4 



-oo 



(2.6) 



in terms of the usual unit polar vectors {e t , e r , eg, e^} in flat space-time. 

By making use of the NP formalism, Teukolsky ( [12], [13]) decoupled the equations 
describing the spin-0, -1/2,-1 and -2 linear perturbations in a general Type D vacuum back- 
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ground. He further showed that, using the Kinnersley tetrad, some of the decoupled equa- 
tions are separable in Boyer-Lindquist co-ordinates in the Kerr background. We will refer 
to the second order ordinary differential equations for the variables r and 9 resulting from 
this separation as the radial and angular Teukolsky equations respectively. The terminol- 
ogy we will use is the following. A bullet • (and a primed bullet •') superscript indicates 
either 'in', 'up', 'down' or 'out' modes, which are standard but, for clarity, are defined in 
Appendix A. Correspondingly, the symbol u* is defined as being equal to u when it is part 
of an expression containing 'in' modes and it is equal to u when the expression contains 'up' 
modes. The separability of the equations describing the spin-1 perturbations implies that 
the electromagnetic potential can be expressed as a Fourier mode sum as 

/+oo +°° +1 

"°° l=\h\m=~l P=±l 

ImwpA'^ = ImjA' + PVlmujA'^ (2.7b) 

where a' m ^p are the coefficients of the Fourier series. The parameter / labels the eigenvalues 
of the angular Teukolsky equation. The parameter P is summed in (2.7a) over the values 
+ 1 and —1, corresponding to two linearly independent polarization states for the potential, 
as we shall see in Section V. The decomposition of the potential into eigenstates of the 
parity operator V is the natural choice because of the invariance of the Kerr metric under 
this operation. 

Chrzanowski [11] was the first author to give analytic expressions, using Teukolsky's 
results, for the linear electromagnetic and gravitational perturbation potentials in the Kerr 
background. He obtained, for the electromagnetic potential modes: 
A m — 

l—m—uiP^fj, — 

= { [1,(6* + 2(3* + O - m;(D + 2e* + p*)} ^R^{r) +1 Z lrmi (9, 0)e"^+ (2.8a) 
+ P [1,(5 + 2(3 + r)-m,(D + 2e + p)] -1^(0-1^(0, 

A U P — 

l-rn-ivPSi, — 

= [p*~ 2 [-n,(5 - 2a* + it*) + m,(A - 2 7 * + p*)] +1 BZ,(r)- 1 Z lmu {9 , 0)e"^+ (2.8b) 

+ Pp*~ 2 [-n,(5* -2a + n)+ m*(A - 2 7 + p)] +1 R^(r) +1 Z lrmJ (9, 0)e— f } 
where 

C_1 \m+l 

h Z lmu {6, 0) = 1 j= h S lmuj (9)e +m t (2.9) 

V 27T 
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The only difference here with Chrzanowski's original expressions is an overall change in the 
sign of m and u>, justified by the sum over m and integration over u. 

It is immediate from the radial and angular Teukolsky equations that their respective 
solutions hRimu)ij) and hSimu}{6) satisfy the following symmetries: 

hRimUr) = &~ h - h R* lmu {r) (2.10a) 
h Rim^(r) = hRlm-^r) (2.10b) 

and 



hSlmuj(0) 



\l+m 



hSl moJ (0) — 

Note the symmetry property 



\h+m 



hSl—m- 



>(0) 



■1) 



l+m 



l—m—w 



A? 



(2.11a) 
(2.11b) 
(2.11c) 

(2.12) 



which is easily obtained from (2.10b), (2.11c) and (2.8). It then follows that the Fourier 
sum coefficients must satisfy the following condition so that the potential remains real: 

<Lp = (-lY^Pal^p (2.13) 

The Maxwell scalars <p' h can be decomposed in a Fourier mode sum analogous to that in 
(2.7a). Following [14] and [10], the Maxwell scalars modes can be expressed as 



-m—uj'T-l 



I— m— 
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9 Din 7 „-iwt 
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lBi muJ up a -iut 



[2.14) 
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(2.15) 
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where iBf muJ = ~i\f muJ +4maa; — 4a 2 u; 2 and is the eigenvalue of the angular Teukolsky 

equation. We are using the definitions of the operators 

CV* =d e ± Q + ncot9 and X>n = <9 r =F + 2n — — — (2.16) 

where Q = —auj sin 9 + raj sin 9 and K = (r 2 + a 2 )cj — am. We use the convention that the 
upper and lower symbols inside braces go with the upper and lower signs in the equation. 

The asymptotic behaviour of the NP Maxwell scalars, separately for outgoing and ingoing 
waves, in the limit r — > +oo is 



(2.17) 



It is therefore the scalar the one with the asymptotically dominant behaviour for the 

upgoing [ingoing] waves. The above asymptotic behaviour (2.17) was originally obtained by 
Newman and Penrose [9] and is commonly referred to as the peeling theorem. 

The Teukolsky-Starobinskii identities relate solutions, radial or angular, of opposite he- 
licity h. These identities are ( [15]): 

DoDo-lRlmu) — -^+lRlmuj £o£l+lSlmv = lBi muJ -lSl mu j (2.18a) 

AP I> A + i^ ma) = 2iBf muJ ^iRi rnu j C Q C x -iSi m u = iBi mw+ iSi mul (2.18b) 

The angular Teukolsky equation has singular points at x = ±1, where x = cos 9. By using 
the Frobenius method it can be found that the solution that is regular at both boundary 
points x = +1 and —1 is given by 

hSirwix) = (1 - x)l m+/l l/ 2 (l + x)^~ h ^ h y lmw (x) (2.19) 

and the function h,yimu>( x ) behaves close to the boundary points as 

oo 

fcW=E4!) (lTx) n for z->±l (2.20) 

n=0 



The solutions hZimui are called the spin- weighted spheroidal harmonics (SWSH). In the 
spherical limit a = 0, the SWSH reduce to the spin-weighted spherical harmonics hXim- 



s 



We give here two properties ( [16], [17] and [18]) satisfied by the spin- weighted spherical 

harmonics which we will need later on: 

i 

-iY lm (0A) + iYL(0A) = (2.21) 

m=— I 

and the "addition theorem" 

1 21+1 

h Y lm (0A)hYL(0'A') = -±-fl(cos 7 ) (2.22) 

m=—l 

where 7 is defined by cos 7 = cos 9 cos 9' + sin^sin^'cos(0 — 0'). 

The classical, electromagnetic stress-energy tensor can be expressed in terms of the NP 
Maxwell scalars as 

v)\ 0+10+1 

1 (2.23) 

- 4000-i^^m^) - A(f) +1 (f)* l {ll m u) + 20 +1 0!l 1 m M m^| + c.c. 

in the case of absence of a charge current. By virtue of the Maxwell field equations, the 
stress-energy tensor (2.23) satisfies the conservation equation T^.^ = 0. 

Note that it follows from (2.5) that all the pairs of null tetrad vectors appearing in the 
different terms in (2.23) remain invariant under the parity operation, except for the ones 
that have a factor containing 0o together with either 0_i or + i, which change sign. That 
is, a pair of null vectors e^e^) appearing in (2.23) with a factor 0^0^/ changes under the 
parity operation as 

V (e {a) e {b) ) = (-l) h+h ' e {a) e {b) (2.24) 

The variable T refers to either the potential components or the NP scalars 0^. The 
Fourier series expansion for either the potential components or the NP Maxwell scalars may 
then be expressed as 

r+00 

r = I du'a' lmujPlmu>P r (2.25) 

ImP 



which may be re-arranged as 



r+oc 

T* = 2_] / du;* (a' muJ pi mL jpT' + (-l) l+m Pa'^p^m^pT') 

ImP Jo 



[2.26) 



We can now use the symmetry relations 



VimuipA' — (—l) l+m i- m -upA°* — PimupA' (2.27a) 



\l+m+l+h 



v lm ^- h = (-ir m+i+ v m ^0r (2.27b) 
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The equations above for T are equally valid for and In particular, we only need to 
apply the operator K£ to an equation for in order to obtain the corresponding equation 
for (ph- However the last step in (2.27a) has no equivalent for <ph in (2.27b). The reason is 
that 

K 'Plnu.A; OC K£ 1-m-uA*; = (2.28) 

as can be checked; that is, this term is pure gauge. Hence the fact that i muJ (f)^ = K£ im^pA'^ 
does not actually depend on P which is why P is not a subindex of the NP scalar modes. 
The potential is real whereas the field components are not, as seen in equations (2.29) and 
(2.30) below. From (2.26) and using (2.27a) and (2.27b) for the potential and the field 
respectively, we have 



r+oo 

^1 = ^2 falTnuPlrnupA'p + P a'^ p V i mu) P A'*) 

lmP J ° 

r+oo 

= / du' (a' mu)P imujpA^ + a'^pim^pA'*) 

lmP J ° 



(2.29) 



and 



/+oo 
dw' {aL.Pim^l + (-l) h+1 PaZ.pVi m ^l*) (2.30) 



III. QUANTIZATION OF THE ELECTROMAGNETIC POTENTIAL/FIELD 

The abundance in the literature of the quantization of the scalar field in a curved back- 
ground is in sharp contrast with the scarce treatment of the quantization of the electro- 
magnetic -or gravitational- field in such a background. CCH did quantize both the elec- 
tromagnetic and the gravitational fields in the Kerr background. They used a canonical 
quantization method, which is the one we have chosen to use in this paper. 

We quantize the field by promoting a' muJ p and a'^p to operators a* muJ p and a'^p respec- 
tively. Cohen and Kegeles [19] and Wald [20] have shown that the theory may be expressed 
in terms of one single NP complex scalar, which represents the two radiative degrees of free- 
dom of the electromagnetic perturbations. If we introduce expansion (2.30) for the NP 
scalars into T 00 given by (2.23), we then obtain a hamiltonian which is a superposition 
of independent harmonic oscillator hamiltonians, one for each mode of the electromagnetic 
field. From the standard quantization of the harmonic oscillator, we know that the operators 
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®imujp an d <hmu>p mus t satisfy the commutation relations: 



a 



\mujPi Q'Vm'w'P 1 
ImuiPi a l'm'u>'P' 



5 {to - u}')5u>5 mm >5pp' 



a lmuiPy a l'm'ui'P' 



(3.1) 



These commutation relations are satisfied provided that the orthonormality conditions 



Im-LopA^, I'm'uj'P'Aa 



S,,'5u>S mm >S(uj - lo')8 



pp> 



ImupA** , I'm'ui'P'Aa ) — 



(3.2) 



are satisfied, where S is any complete Cauchy hypersurface for the outer region of the 
space-time and where the Klein-Gordon inner product is taken as 



{lJ>a,<Pa)i 



i J d'w (v a * v^ Q - <^*v M c + ^v Q v Q * - r^ a ^ a ) (3.3) 



The inner product (3.3) has the same form as the one taken by CCH. However, CCH give an 
expression for the stress-energy tensor which includes a factor 4ir in (2.23), corresponding 
to unrationalized units. If unrationalized units are used, then a factor 4tt should also be 
included in the inner product (3.3). 

Constants of normalization are to be included in front of the radial functions so that 
the potential modes (2.8) satisfy the orthonormality conditions (3.2) given the asymptotic 
behaviour of the radial functions in (A4). We find that the constants of normalization are 
given by: 

1 



in 1 2 



1^-1 



\K p i\ 2 



2 5 cu 3 



7T 



2%|9I| 2 cu(r2 +a?) 



up 1 2 



\NZ | 2 



up,mc 



j-)Up,mc 
+l n imu) 



uj{r 2 + + a 2 ) 



where 91 = iK + + (r + — r_)/2 and K + = K(r + ). We have chosen _\R 



Imu) 



(3.4a) 
(3.4b) 

(3.4c) 

1 and 



up,inc 
Imui 



1 in equations (3.4a) and (3.4b) respectively. The constant of normalization \N^\ 



is calculated as indicated with the use of (A10). It is therefore the constant of normalization 
that corresponds to using the radial function (A4b) when setting -i-R^^ nc = 1, which is the 
actual normalization we have used in the numerical calculation of the 'up' solutions. The 
NP scalars are therefore assumed to include the constants of normalization (3.4). That is, 
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the NP scalar modes z maJ 0* are to be calculated from expressions (2.14) and (2.15) with 
the inclusion of the appropriate constant of normalization (3.4), while the radial functions 
remain unaltered. 

It can be checked that the set of modes {imupA™, imuopA^} forms a complete set of or- 
thonormal solutions to the Maxwell equations in the outer region of the Kerr space-time. 
Similarly, it can be checked that {; mw pA° ut , imupA^ 11 } also form a complete set. We may ex- 
pand the electromagnetic potential by using the complete set of solutions {imupA™, i mu jpA np } 
and then quantize it as: 

" + 00 



Ap, = / da; [a^pimupA™ + a^ w pi mu)P A^* J + 
imp Jo V 7 

r+oo . . 

+ ^2 ^ ( ^IrLp WP^f + Kmu,PlmuipA^ P * ) 

;™d^0 V 7 



(3.5) 



IraP 

Alternatively, we could proceed exactly in the same manner but using the complete set 
of solutions {imuipA ^, i m(J pA^ wn } instead. The result is then: 

Ap = y2 ^ ( ^mujPlrnuipA ^ + a^ptapA™'* ) + 

; d JO X 7 

r+ oo , . (3-6) 

i \ ^ / A,~, I f, down A down , ~downf /| down* \ 

{^ImuPlmuPAp + O-lmujPlmujpA^ J 

ImP ^° 

The asymptotic behaviour in terms of the advanced v and retarded u time co-ordinates 
of the electromagnetic potential and NP scalars for the 'in' and 'up' modes is the same as 
the one exhibited by the modes (A2). The same applies to the asymptotic behaviour of 
the 'out' and 'down' modes exhibited in (A16). Accordingly, the operators a^p, a^ p , 
®°mlp an d ofmI,p are creation operators of particles incident from past null infinity X~, past 
horizon 7i~ , future null infinity X + and future horizon 7i + respectively. 

Since the 'in' and 'out' modes are only defined for uj non-negative, they have non-negative 
energy as measured by an observer following the integral curve of ky, by virtue of (A3). 
Similarly, the 'up' and 'down' modes, defined for Co non-negative, have non-negative energy 
with respect to observers following the integral curve of k^. 

We may now construct the stress-energy tensor operator from either the potential op- 
erator (3.5) or (3.6). It is well-known that the stress-energy tensor as an operator does 
not have a well-defined meaning. It suffers from ultra-violet divergences and its expecta- 
tion value when the field is in a certain state must be renormalized. There are several 



12 



techniques for renormalization. The point-splitting technique consists in starting from each 
quadratic term in the stress-energy tensor and temporarily displacing one field point, thus 
forming the object ^T Q/ g(x, x r )J , which is finite. Specific divergent terms, gathered in the 
bitensor T^(x, x'), which are purely geometric and thus independent of the quantum state, 
are then subtracted from (r a p(x, x')\ . The end result is obtained by finally bringing the 
separated points together: 

(f Q/3 (x))* n = Km (jf a p(x,x')y - T^(x,x f )j (3.7) 

It is this renormalized expectation value of the stress-energy tensor (RSET) that is the source 
in Einstein's field equations in the semiclassical theory. Christensen [21] has explicitly 
calculated the divergent terms T^ v by using covariant geodesic point separation. Jensen, 
McLaughlin and Ottewill [22] calculated a linearly divergent term for the spin-1 case, which 
was not explicitly given by Christensen. The reason being that this term does not have to 
be included when an average is taken over the covariant derivative of the biscalar of geodetic 
interval <r M and — cx M , as performed by Christensen. 

Before we start a description of the various physical states of the field, we give an im- 
portant result found by Unruh [23] and further established by [24] and [25]. The result is 
that a 'particle detector' will react to states of the field which have positive frequency with 
respect to the detector's proper time. If a certain observer A makes measurements relative 
to a certain vacuum state |S), then he or she measures a stress tensor 

when the field is in a certain state 1^). 

IV. VACUA STATES 

The properties of the various states described in this section have been obtained in the 
literature for the scalar case, except where explicitly indicated otherwise. 

The Boulware vacuum state, denoted by \B), is defined in Schwarzschild space-time as 
the vacuum that corresponds to quantizing the field with normal modes that have all positive 
frequency with respect to the space-time's hyper surf ace-orthogonal timelike killing vector 
ky. This state respects the isometries of Schwarzschild space-time. Since it is the static 
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observers SO the ones that move along integral curves of k j?, from Unruh's result stated in 
the previous section it follows that these observers will make measurements relative to the 
Boulware vacuum \B). Candelas [6], based on conjectures made previously by Christensen 
and Fulling [26] , has found that the RSET when the scalar field is in the Boulware vacuum 
is zero at both X~ and X + in the Schwarzschild space-time. Candelas also found that 
the RSET, close to the horizon, when the field is in the Boulware vacuum diverges and 
corresponds to the absence from the vacuum of black-body radiation at the black hole 
temperature appropriately red-shifted. The Boulware vacuum is therefore irregular at lir 
and 7i + . 

In Schwarzschild space-time the Boulware vacuum may be associated with the field ex- 
pansion in terms of either the 'in' and 'up' modes or the 'out' and 'down', both pairs of sets 
of complete modes defining the same vacuum \B). We can perform a similar expansion for 
the electromagnetic potential in the Kerr space-time. The past Boulware state is defined by 

I / 

B-) = 



- up 
a lmuiP 



corresponding to an absence of particles at 7i and X . We can also define the future 
Boulware state, as that state which is empty at X + and 7i + : 

aZip\B + ) = 
B + ) = 



-down 
l lrruuP 



The Bogolubov transformation between the pair of operators ctf^p and a^p and the pair 
af^p and af^p is non-trivial: the expression for a^p 



*upf 
a lmujP 



in terms of 'out' and 'down' 



, , • -downt 

operators contains a t _ m _ u P 



- outf 

a L-m,-ui,P 



for modes in the superradiant regime. This implies 
that the past Boulware state contains both outgoing and downgoing superradiant particles, 
and is therefore not empty at X + and 7i + . This flux of particles out to X + corresponds 
to the Starobinskii-Unruh effect. Similarly, the future Boulware state contains ingoing and 
upgoing superradiant particles, and is therefore not empty at X~ and 7i~ . As ajjj^p, when 
expressed in terms of 'out' and 'down' operators, contains the creator operator a°_ m _ UJ P , it 
is not possible to construct a state which is empty at both X~ and X + , unlike the situation 
in the Schwarzschild space-time. 

From the definitions (4.1) and (4.2) together with the relations (A17), the past and 
future Boulware states are obtainable one from the other under the transformation (t, <ft) — > 
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(—t, —0). Because the two states are not equivalent, it follows that neither is invariant under 
this symmetry of the Kerr space-time. 

The defining features of a Hartle-Hawking state ( [27]) is that it possesses the symmetries 
of the space-time and that it is regular everywhere, including on both the past and the future 
event horizons. Kay and Wald [28] have proven that for any globally hyperbolic space-time 
which has a Killing field with a bifurcate Killing horizon there can be at most one state 
with the above features. Kay and Wald have further shown that for the Kerr space-time 
this state does not exist. Rindler and Schwarzschild space-times are covered by Kay and 
Wald's theorem; in Rindler space-time this state is clearly the Minkowski vacuum. 

In Schwarzschild space-time the state \H) corresponds to quantizing the field with up- 
going normal modes which on 1~C~ have positive frequency with respect to the Kruskal 
co-ordinate U = — e~ K+u and with ingoing normal modes which on 7i + have positive fre- 
quency with respect to the Kruskal co-ordinate V = e K+v . Candelas [6] showed that the 
state defined in this manner is regular on both the past and future horizons. He also found 
that the RSET at infinity when the field is in the \H) state corresponds to that of a bath 
of black body radiation at the black hole temperature Th = k+/(2it). The Hartle-Hawking 
state models a black hole in (unstable) thermal equilibrium with an infinite distribution at 
the Hawking temperature. 

I- \ H ~ B 

From the above results and from the previous section we know that ( T a p j is thermal 

both for r — > r + and for r — > +oo. Christensen and Fulling conjectured that this is the 

case everywhere. However, Jensen, McLaughlin and Ottewill [29] numerically showed that 
I- \ H ~ B 

( T a p j deviates from isotropic, thermal form as one moves away from the horizon. 

Candelas showed that the RSET close to the horizon when the field is in the \B) state 
diverges like minus the stress tensor of black body radiation at the black hole temperature, 
and that it must tend to — (r a p*j , due to (3.8) and to the regularity of \H). Analogously, 
Unruh [23] showed that in flat space-time and when the field is in the Minkowski vacuum, a 
Rindler observer RO will also see a bath of black body radiation at the Hawking temperature 
of a black hole with surface gravity k+ = aa, where a is the RO's acceleration and a is the 
lapse function in Rindler space. 

Frolov and Thorne [30] defined a new "Hartle-Hawking" state \FT) invariant under the 
symmetries of the Kerr space-time by using a variant of the t] formalism (which employs 
non-standard commutation relations for the creation and annihilation operators). They 
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proved that the RSET when the field is in the \FT) state is finite at the horizon but 
that, at least for arbitrarily slow rotation, it is equal to the stress tensor of a thermal 
distribution at the Hawking temperature rigidly rotating with the horizon. This suggests 
that it becomes irregular wherever h,%> is not timelike, that is on and outside the speed- 
of-light surface. Ottewill and Winstanley [31], however, proved that although \FT) has 
a Feynman propagator with the correct properties for regularity on the horizons, its two- 
point function is actually pathological almost everywhere, not just outside the speed-of-light 
surface. Only at the axis of symmetry, where all the modes in the two-point function for 
the scalar field are evaluated for ui = u (i.e., m = 0), it does not suffer from this pathology. 

Frolov and Thorne claim that close to the horizon ZAMOs make measurements relative to 
an unspecified Boulware vacuum. They also claim that, when the field is in the state \FT), 
ZAMOs measure close to the horizon a stress tensor equal to that of a thermal distribution 
at the Hawking temperature rigidly rotating with the horizon. 

Duffy [32] modified the Kerr space-time by introducing a mirror and constructed a state 
\Hm) for the scalar field that is invariant under the isometries of the modified space-time. 
He then showed that \Hm) is regular everywhere in the modified space-time if, and only 
if, the mirror removes the region outside the speed-of-light surface. He constructed another 
state, |-B.m), invariant under the isometries of the modified space-time and empty on both 
the past and future horizons. This is the state that RROs make measurements relative 
to in the modified space-time. He also numerically showed that when the field is in the 
\Hm ) state the stress tensor measured by a RRO is, close to the horizon, that of a thermal 
distribution at the Hawking temperature rigidly rotating with the horizon. 

CCH defined a new Hartle-Hawking-type state, which we will hereafter denote by 
\CCH~). This state is obtained by thermalizing the 'in' and 'up' modes with respect 
to their natural energy. Ottewill and Winstanley [31] showed that this state is, however, 
not invariant under the symmetry transformation (t, <p) — > (— t, — 0) of the space-time. They 
further argued that the RSET when the scalar field is in the \CCH~) state is regular on 
the future horizon but irregular on the past horizon. We must note that these results were 
derived in [31] based on a stress-energy tensor for which the t8- and ^-components are 
identically zero. We shall see in Section VI that although this is indeed the case for the scalar 
field, which is the case they considered, this is most probably not true for the electromagnetic 
field. A similar state could be constructed by applying the transformation (t, <fi) — > (— t, —(f)) 
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to the state \CCH ). This state, suitably named \CCH + ), would then be irregular on the 
future horizon and regular on the past horizon. In Section VII we will investigate the form 

/ „ xCCH--B- 

close to the horizon of ( ) for electromagnetism. 

\ / rcn 

Finally, Unruh [23] constructed a state in the Schwarzschild space-time by expanding the 
scalar field in modes that are positive frequency with respect to the proper time t of inertial 
observers in X~ and modes that are positive frequency with respect to the proper time of 
inertial observers close to IHr . It is possible to construct a state in the Kerr space-time 
with the same positive-frequency mode definitions. We call this state the past Unruh state, 

\u-). 



V. EXPECTATION VALUE OF THE STRESS TENSOR 



A. Analytic expressions 

In this subsection we derive new, analytic expressions for the expectation value of the 
stress-energy tensor when the electromagnetic field is in quantum states of interest in the 
Kerr space-time. It is well-known that there is an operator-ordering ambiguity in the tran- 
sition from classical to quantum theory: the classical term (fih^h' should be quantized to the 
symmetrized form ^</>h</>jj/ + (fih'&hj We will calculate separately the expectation value 
of the two quadratic terms in the symmetrized form. 

It is straight-forward to check that in the past Boulware state we have 



(r+oo p- 
Jo Jo 

= {-l) h+h 'v((B-\^ h \B-)\ 



(5-1) 



The sign (— \) h+h ' is precisely the same sign appearing in (2.24). This implies that if 
the quadratic terms in the expression (2.23) are quantum-mechanically symmetrized when 
promoting the NP scalars to operators, then the expectation value in the state \B~) of the 
stress-energy tensor will be invariant under parity. 

In order to calculate the expectation value of the quadratic terms 4>h<t>\i and (f) h i4>h in the 
past Unruh state we are going to make use of the expression calculated in [30] which gives 
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the past Unruh state in terms of the past Boulware state: 



\u~) = n ^ P exp( e -^a^ P a£; p ) \b- 



(5.2) 



ImCbP 



where Ci muJ p are normalization constants and a"^p are creation operators in the left hand 
region of the extended Kerr space-time. We will also make use of the following expression 
in [33]: 



S(r, 



(coshr) 1 exp ^— a+aLe 2l<?i tanhr J exp ^— (a+a + + aLa_) ln(coshr)j exp (a + a_e 2i *tanhr) 



where S^r, </>) is the two-mode squeeze operator 



S(r, 0) = exp ( r(a + a_e" 2 ^ - a ] + a}_e 2i,f> ) 



(5.3) 



(5.4) 



and the independent operators a + and a_ satisfy the standard commutation relations. By 
using (5.3) and (5.4) we can re-express (5.2) as 



JmP 



U ) = exp < ^2 / du; In C imw p + ln(coshr^) >e A |fi ) 



(5.5) 



with 



up'f ~up *up' 
ImujP ~ a lmujP a lmuiP 



and r,~, 



tanh" 1 (e"^ /K +) . 

(5.6) 



Since (e A ) = e A , the normalization (U \U ) = 1 impli 



ICS 



JmP 



exp 1 / dcJ ln QmuP + hi Cfma,p + 2 ln(cosh r<; 



Using now the Baker- Campbell-Hausdorff equation ( [34]) we find that 



(5.7) 



A * up 
e a lmuiP e 



HmwP 



coshr^ + aj^jp sinhr £ 



(5.* 



and finally 



( U | O'lmuiP^lKn'uj'P' \U ) — o 



. 7T0>\ 

coth ( — - 1 



PP , 



( ^ I «lS,p a £'o,'P' \ U ) = = {U I a^pa^p, |c/ > 



(5.9a) 
(5.9b) 
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It immediately follows that 

{u-\m\,\u-) = 



■y poo 



ImP 



-l) h+h V(i muJ cp h IrmjVhi 



i up* 



up\ 



+ 

POO 

+ dcu 
Jo 



Imuj'Vh Imuj'Vh' 



i up* 



<up\ 



Imuj'rh IrwYy 



, ( iru\ 

cotb {-J 



+ 



(5.10) 



Note the minus sign in the second term in (5.10). Its presence may seem a bit surprising at 
first but, as we shall now see, it is precisely this sign that causes the expectation value of the 
stress-energy tensor in the past Unruh state to adopt a more familiar form by having all 'up' 
terms multiplied by a coth factor. This is already clear from looking at (5.10) and realizing 
that when quantum-symmetrizing the classical expression 4^4^* the terms without a coth 
factor will cancel out. 

The following identities are therefore immediately satisfied 

B~ 















[k, A{ 



E 

ImP 



lmui ( t ) hlmu)4 > h' 



-1 



V 



(5.11) 



and hence 



U--B- 



(5.12) 



as it should be. 

When the classical term 



^h^h' + c.c.) is quantized to the symmetrized term 
h<p\i + 4>\'4>h) /2 + h.c, where the symbol '/i.e.' stands for hermitian conjugate, we obtain 
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the following real, parity-invariant expressions in the past Boulware and past Unruh states: 



B- 




\ lmLO<P h i + (,-lJ y\lmw(p h lmw<Ph' ) 



lmu)<Ph 

U~ 



+ 



(5.13a) 



c.c. 



/up iup 
Imoj'Ph lmu)<Ph' 



p* + (_l)M-*>( Jmw ^ Imw ^)] coth (™) 



(5.13b) 



c.c. 



ImoiClmw 1 ?!'™ + ( — l) h+h Vdmufthlmuftfr* 

Note that in the above expressions we have been able to complex conjugate the mode 
functions that are operated on by V because of the existence of the +c.c. terms. This 
immediately leads to the following real, parity-invariant expressions for the stress-energy 
tensor in the past Boulware and past Unruh states: 



B T^ u B 



>*] + (-irp (t,a 



1 ImP ^ 

4e 



rr\) }+ 

in*] \ 
h \) 



(5.14a) 



ImP 



dcu coth {^\ [T, u [ lmu ^, lm ^r\ + (-1)^ (T, v [ lm ^f, i^T\) }+ ( 5 - 14b ) 



+ J duo [t^ [ lmu </%, lmu <FF] + (-ifv (t^ [ 



im nil* 



]) 



We use the obvious notation that T^ u [i mLU (f>* , imu^'C} denotes the general expression for the 
stress-energy tensor (2.23) where the scalars <ph have been replaced by the modes imu^h ■ We 
will also use the symbol irn^T* u to refer to T^ v [imu<f>%, imw^V]- The variable $ is defined so 
that (— 1)® is equal to -1 if one index of the component of the stress tensor is 9 and the other 
one is not, and it is equal to +1 otherwise. Note that the sign (— 1)^ appears in the above 
expressions instead of (— \) h+h ' by virtue of the change under the parity operation of the 
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coefficients of the quadratic field operators that appear in the expression for the stress-energy 
tensor, as seen in (2.24). 

To our knowledge, these important expressions for the electromagnetic field in the Kerr 
space-time have only been given so far by CCH. Their expressions, however, differ from ours 
in that they do not include the terms explicitly containing the operator V. As a result, the 
expressions for the expectation value of the stress tensor that CCH give are not symmetric 
under the parity operation V : (6, <p) — > (tt — 6, <p + tt), even though this is a symmetry of 
both the Kerr metric and the Maxwell field equations. 

We initially used CCH's expressions in our numerical calculations and the numerical 
results we obtained were clearly not symmetric under the parity operation. We also showed 
analytically that at least for one particular instance the RSET is not symmetric under V. 
We outline here the derivation. 

From equation (2.20) and the two equivalent expressions for i mul 4>o given in (2.15) it 
follows that |/ m a;0o| 2 is only non-zero at 9 = 0, ir if m = and its value, for the 'up' modes, 
in that case is 

\i,m=o^(pl P (r,9 = 0)| 2 - \i, m =Q, w 4>l v {r,6 = tt)| 2 = 

-4ia\Nl\\ 2 W[ +1 R, -iR*]^ =Q J-ian=o,i,m=o,u, +ian=cM,m=o,J ( 5 - 15 ) 

By virtue of the property (Al9a), the corresponding value for the 'in' modes is obtained by 
merely changing the sign. 

Since T^™ is a purely geometrical object and the metric is invariant under V, this divergent 
stress tensor must also be invariant under V. From this property together with equations 
(2.23) and (5.15) we find that CCH's expressions for the expectation value of the stress 
tensor yield: 

f ee (r,e = Q)) U -(f 0e (r,6 = 7r)) U = 

ren 

(5.16) 



I rcn 

oo 

, / TIUI \ 

X 



coth (_) - 1 



< 2 3 za|A r ^| 2 iy[ + i-R, -i-R*]^ =0i ^ +i a n=o,;,m=o^ i -i\m=o,w + 2acu 
(r 2 + a 2 )\Bl m=0 uj y -iAi )TO=0)W - 2ouj 

where we have made use of the Teukolsky-Starobinskh identities (2.18) to relate the SWSH 
for h — +1 with the one for h — — 1 in the limit x — > +1. 
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It is easy to check that iW[ + iR, -iR*]^ p 



m=0,ui 



l, m=0, u 



up,tra 
m=0,u; | 



juj > as long as 



uj > 0, and therefore the integrand in (5.16) is non-negative for uj > 0. It follows that 
(Tqq\ is not symmetric under V at the axis when CCH's expressions are used. 

\ / rcn 

On the Schwarzschild background, CCH's expressions for the expectation value of the 
electromagnetic field are straight-forwardly invariant under parity: when a = 0, the angular 
mode function does not depend on u> and the radial mode function does not depend on m. 
Therefore, applying the transformation {9 — ► n — 0) on any mode in their expressions is 
equivalent to performing (m — > —m) only, by virtue of the relations (2.10b) and (2.11b) 
and the reality of the stress tensor. The presence of the sum Y^ m =-i m their expressions 
guarantees their invar iance under V. Such a straight-forward reasoning does not follow in 
the Kerr background because applying the parity operation on CCH's expressions implies a 
change in the sign of uj as well as in the sign of m. Whereas the sum over m is symmetric 
with respect to m — 0, the integration over the frequency is not symmetric with respect to 
uj = 0. 

CCH's expressions are, however, valid when applied to the scalar field case in the Kerr 
background. By comparing the scalar field's expectation values and CCH's expectation 
values with the symmetrized versions (5.14) for the states \B~~) and \U~), we can give 
analogous symmetrized versions for the states \FT) (not included in CCH as this state was 
only defined later) and \CCH~): 



FT 



T, 



FT 



2 ^ 



ImP 



du coth 



7TUJ 



(5.17a) 



+ / du; coth ( 



{CCH' 



T 



/ HUJ 

CCH 



+ 



do) coth 

o 

du; coth I 



7TUJ 
K + 
1XUJ 



2 ^ 

T^ v [imcj Imw C ] 



CD 



ImP 



i up /A U P* 1 



C*D 



(5.17b) 



CD 
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B. Polarization 

In this subsection we will give a physical interpretation of the non-parity term and the 
parity term appearing in the expressions (5.14) and (5.17) for the expectation value of the 
stress-energy tensor in different states. We denote by parity term in a certain expression a 
term that explicitly contains the parity operator V, and by non-parity term one in the same 
expression that does not explicitly contain this operator. 

It is clear from the classical expression (2.8b) for the 'upgoing gauge' potential i m u>A upiJl 
that this potential only contains the null vectors n and m. The parity term PVi mu} A np11 
in (2.7b), because of the transformations (2.5) of the null base under parity, contains the 
vectors n and m*. The potential i m uipA up ^ therefore contains the vectors n, m and m*. 
However, it is well-known ( [35]) that only two of them are physically significant at radial 
infinity, where the space-time is flat. 

It is in the limit for large r that the physical meaning of the various vectors becomes 
clear. We know that in flat space-time an electric field mode of positive frequency with 
radial and timelike exponential behaviour e-M*- r ) (which we shall see is the asymptotic 
behaviour in our case) that is proportional to the vector (eg + ie^) [(eg — ie^)} possesses a 
negative [positive] angular momentum and we thus say that it is negatively [positively] polar- 
ized. If the mode is instead of negative frequency, the sign of the angular momentum changes 
and then an electric field mode with radial and timelike exponential behaviour eT^ t ~ r > pro- 
portional to (eg + ie$) [(eg — ie^)} is said to be positively [negatively] polarized. Therefore, 
according to (2.6), an electric and a magnetic field modes of positive frequency that are 
proportional to the vector m[m*] correspond, in the flat space limit, to a negative [positive] 
polarization, whereas the vectors I and n correspond both to neutral polarization. 

We know from (2.28) that the parity term is pure gauge and therefore the contribution 
to the NP scalars from the term with m* in the potential is zero. Only the terms with n 
and m in the 'upgoing' potential contribute to the NP scalars. In particular, it is immediate 
from expressions (2.4) for the Maxwell scalars that only the term in the potential im^pA^ 
that contains the vector m contributes to i m u<l>+x whereas only the term with n contributes 
to imui^i- Both, terms with n and terms with m, contribute to / ma ,0o P . This implies that 
the positive-frequency modes imu^i, imw^i an d imu^ are obtained from terms in the 
potential that, in the flat space limit, are neutrally-, negatively- and both neutrally- and 
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negatively- polarized respectively. The parity term in the potential in (2.7b), which is of the 
opposite polarization to that of the non-parity term, does not contribute to any NP scalar 
mode because it is pure gauge. The opposite polarization to that of the non-parity term in 
the potential contributes to the NP scalars through the negative-frequency modes when the 
integration is over all frequencies, as in (2.25). In the expression (2.29) for the potential or 
(2.30) for the NP scalars, in which we have rid of the negative-frequency modes, the opposite 
polarization appears via the complex-conjugate term or the parity-term respectively. 

Even though all three Maxwell scalars appear in the classical expression for the electro- 
magnetic stress tensor, due to their different asymptotic behaviour (2.17) for large r, the 
terms in the stress tensor (2.23) with i mU )(t^i predominate in this limit. That is, the radiation 
field components of the stress tensor are calculated in the flat space limit from modes in 
the potential (2.25) which for positive [negative] frequency correspond to a negative [positive] 
polarization. Note that the complex-conjugation of NP scalars in the stress tensor does not 
change the polarization of the field since it is merely a consequence of the fact that the null 
tetrad contains complex vectors, and does not imply the complex-conjugation of the tensor 
field components F^ u . 

We give here expressions for the 'upgoing' field modes in the limit for large r, where 
the space-time is flat. We wish, however, to obtain expressions for the fields that are real 
mode by mode. We will therefore not calculate them from the potential modes im^pA^ in 
(2.8), since the transformation (m, u) — > (— m, —uj) has been applied to the parity term in 
(2.7b). Even though the potential A^ is obviously real, the potential modes imuipA^ are 
not. Instead, we will calculate field modes from the potential modes: im^A^ + i mu) A np * . The 
result is: 

/ > 2 1 N np I 
pup _^ V l iV +ll v 
ImcuP-n — ► 7= X 

V2r 

x |_iYj_ m _ w +iKi_ m __ UJ e (e e + ze J + _iY,_ m _ w +iK l _ m _ UJ e (eg - le.^ 



-^ 2 z|iV up 

ImujP-D — > 



x [-lYi-m-Lu +\ri l _ m _ UJ e ' (e e + xe^) - -i^^^ +iK l _ m _ UJ e (e e - ie q 



-oo 



(5.18) 

where we have used the fact that in flat space we can replace hZimu) by bXimu- It is worth 
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noting that even though the second term in either (z ma; A^ p + im^A^*) or (2.7b) does not 
contribute to the NP scalars it does contribute to the fields. The reason is that the NP 
scalars must be calculated from real fields. It does not make physical sense to consider the 
contribution to the NP scalars from a non-real field, such as the second term in the above 
expressions for the field. We say that this term is 'pure gauge' in the sense that it does not 
contribute to the NP scalars even if it does contribute to the physical fields so as to make 
them real. The large-r asymptotics for the NP Maxwell scalars in terms of the electric and 
magnetic fields are easily obtained: 

0-i -> — 7= (E + Bi) (eg + ie^) (r -> +oo) 
V2 

0o -> - (E + Bi) e r (r^+oo) (5.19) 

0+i -> -^j= (E + Bi) (eg - ie<j,) (r -> +oo) 

It is clear that the parity and non-parity terms correspond to opposite polarizations for the 
fields (5.18) (and also for the potential). It is also clear that the only contribution to 
from the 'upgoing' electric and magnetic fields comes from the non-parity term. To leading 
order in r for the electric and magnetic fields both and </>q P vanish, in agreement with 
(2.17). To next order in r, expressions (5.18) and (5.19) must be calculated to include lower 
order terms and it is therefore not valid to conclude from them that the only contribution 
to 0" p and 0q P comes from positively- and neutrally- polarized terms respectively. We have 
indeed seen in the beginning of this subsection that this is not the case. 

The reasoning used so far for the 'upgoing gauge' potential can be applied in the same 
manner to the 'ingoing gauge' potential ; mw v4™. In this case, the potential contains one 
term with the vector I, which is the only one that contributes to imu^+n an d one term 
with the vector m*, which is the only one that contributes to i muJ <f)™ 1 . Like in the 'upgoing 
gauge' case, both terms contribute to z mw 0o P , and the parity term (containing I and m) 
does not contribute to any of the NP scalars. The scalar imuft-i is the one that diminishes 
more slowly in the limit for large r. The radiative components of the classical stress tensor 
is thus calculated in the flat space limit from modes in the potential (2.25) that for 
positive [negative] frequency correspond to positive [negative] polarization. 

So far in this subsection we have looked at the physical meaning of the different terms 
in classical expressions only. We are now in a position to understand the physical meaning 
of the terms in the quantum field theory expressions. The positive frequency modes in 
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(2.25) correspond to the non-parity term in the expression (2.30) for the NP scalar and, 
ultimately, give rise to the non-parity term in the expectation value of the stress tensor 
(5.14). Similarly, the negative frequency modes in (2.25) give rise to the parity term in the 
NP scalars and the parity term in the expectation value of the stress tensor. We therefore 
reach the conclusion that the non-parity terms in expressions (5.14) for the expectation value 
of the stress tensor correspond in the flat space limit to one specific polarization (negative 
in the 'up' case and positive in the 'in' case) and that the corresponding parity terms in 
the same expressions correspond to the opposite polarization. Both the contribution from 
the positive-polarization terms and from the negative-polarization terms are separately real, 
as it should be. We also know, from the end of Subsection VA, that in the spherically- 
symmetrical case a = 0, the contribution to the expectation value of the stress tensor from 
the positive-polarization terms is identical to the one from the negative-polarization terms, 
as one would expect. 

The notable exception to this picture are the 'up' superradiant modes. Indeed, these 
modes have a sign of u opposite to the non-superradiant modes in the same term in the 
expectation value, whether the parity term or the non-parity term. The polarization of the 
'up' superradiant modes is therefore the opposite to the non-superradiant modes in the same 
term, that is, it is positive if part of the non-parity term and negative if part of the parity 
term. Note, however, that the 'in' superradiant modes have the same sign of uj (positive), 
and therefore the same polarization, as the non-superradiant modes in the same term in the 
expectation value. 

When CCH only include non-parity terms in their expressions for the expectation value 
of the stress tensor they are only including one polarization and leaving out the other one 
for the 'in' modes. For the 'up' modes, they are only including one polarization for the non- 
superradiant modes and the opposite polarization for the superradiant modes. In particular, 
when subtracting the expectation value of the stress tensor in the past Boulware state from 
the one in the past Unruh state, only 'up' modes are needed. Neglecting the parity terms is 
in this case equivalent to neglecting positive polarization non-superradiant modes as well as 
negative polarization superradiant modes. That is the case in the calculation of \T^v) 

\ / ren 

close to the horizon in Section VII but, as explained in that section, in this limit the non- 
parity and the parity terms coincide. 

It is interesting to group the terms with the same polarization in the expectation value of 
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the stress energy tensor. Of course, in the case of the difference between the states \CCH~~ ) 
and \U~), which only has contribution from the 'in' modes, the sum of the positive po- 
larization terms coincides with the direct evaluation of CCH's expressions. The negative 
polarization contribution can be obtained by applying the transformation x — > —x. We 
include the plots of the tensor components corresponding to the fluxes of energy and an- 
gular momentum from the positive polarization terms in Figures 1-2. The evaluation of 
the positive polarization contribution to the difference in the expectation value of the stress 
energy tensor between the states \U~) and \B~) requires carefully adding the contribution 
of the superradiant modes to the appropriate polarization. We calculated the positive po- 
larization contribution to these differences of expectation values and plot them in Figures 
3-4. The corresponding negative polarization contribution is, again, obtained by applying 
the transformation x — > —x. The interest of these graphs lies in the region far from the 
horizon. Close to the horizon the irregularity of the state \B~) dominates and thermality 
guarantees symmetry with respect to the equator. 
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FIG. 1: Positive polarization terms of x-A \ Tt r ) 
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FIG. 3: Positive polarization terms of 4- (Tt r ) 




VI. LUMINOSITY 

Let E^ nc ^ and E^ o{ ^ denote, respectively, the energy incident and the energy reflected by 
the black hole at infinity Let _E'( tra ) denote the energy going down across the event-horizon 
of the black hole. The reflection coefficient M.i miU and the transmission coefficient Ti mu) of an 
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incoming wave mode are then denned as the following flux ratios: 

The wronskian relations for the solutions of the radial Teukolsky equation correspond to 
the conservation of energy law that equates the net flux of energy coming in from infinity 
to the net flux of energy going down into the black hole: 

1 - Rjmu; = ^Imcu = 1 - \^Wnu:\ = W[ + lR, -l-R*]^ (6.2) 

where we have made use of the relations (A10) and (A5). For superradiant wave modes the 
fractional gain or loss of energy of an incoming wave mode, i.e., M; mw = l^lmJ 2 , * s greater 
than one while the transmission coefficient Ti mu) is negative. 

The conservation equations V^T^ = can alternatively be written [36] as 

d u (W=S = \V=9 (d,g af s) T a ? (6.3) 

Assuming that the stress-energy tensor is independent of t and 0, like the Kerr metric, the 
fi = t and /i = components of equations (6.3) can be integrated over r to yield: 

T tr = ^~ -r^de (sin 6 f dr'T te ) (6.4a) 



A A sin 9 



r 



^HT-^Le^jJ^) (6 ' 4b » 

where K(9) and L(9) are arbitrary functions. The luminosity when the field is in the state 
| \P ) is given by 



dM 



dt 



dn(f tr .y (6.5) 

/ ren 



where the surface S can be any surface of constant t and r. In the forthcoming the sub index 
A refers to either t or and the subindex X to either r or 6. 

We will now compare some spin-1 results with the corresponding spin-0 results. The 
following results for spin-0 are proven in part by Frolov and Thorne [30] and extended by 
Ottewill and Winstanley [31]: 

f Ae y =0, A/fS =K(6), A/f^Y =L(6) for s = (6.6) 

/ ren \ / ren \ / ren 

where |^) is any state among \B ± ), \U~), \CCH~) or | FT). The last two equations 
in (6.6) are a direct consequence of the first one and of (6.4). It may indeed be calculated 
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directly from the expression for the spin-0 stress-energy tensor that all the radial dependence 
of A imuJ T tr can be expressed as a radial wronskian. It can also be checked that i mLU T™ = 
—imJ^tr for spin-0 so that the only contribution to the luminosity in the past Boulware 
vacuum comes from the superradiant modes: 

\^ \^ / A,., A. T^ U P 



rcn 

1=1 m=l 



A(T tr ) = -2 2^ 2^ / duA^T^ for s = (6.7) 







It is immediately apparent that for the spin-1 case the task of proving analytically whether 
AT tr is constant in r or not is not as straight-forward as for spin-0. When using expressions 
(5.14) for the expectation value of the stress-energy tensor for a spin-1 field, we are interested 
in the following calculation: 

A {imuiTf-jF + 'PlmuiTfJ') = 

+ asin# (-lSimudo-iSimu — +iSi mLU do + ±Si muJ ) \ 

where we have made use of the symmetry (2.11a) in order to relate terms that contain 
i-iSfmu, + v -i S LJ to terms that contain (+i£f maj + V +1 Sf m J. The corresponding result 
for the 'in' modes is equal to (6.8) with a change of sign, by virtue of (6.2) and the property 
(Al9a). We can get rid of the derivatives in (6.8) by using the Teukolsky-Starobinskii 
identities and the angular Teukolsky equation. But even when doing that, it is not possible 
to express the last term in (6.8) in terms of -\Sf mul and + iSf muJ only. As a matter of fact, 
when evaluated at the axis of symmetry 9 = or it, only one term containing _i5*; 2 mw and one 
term containing +iSf muJ appear, neither of which is constant in r. It is therefore apparent 
that if we wish to prove that A (i mu jT^F + Vi miJ jT^F) is constant in r, or otherwise, we must 
then somehow relate -iSf mu} to +iSf mu} . It follows from the symmetries (2.11) that we can 
only relate at the same point a term whose spherical dependence is _iSf muJ to a term whose 
spherical dependence is + iSf mLU by applying the transformation (m,cu) —> (— m, —uj) to one 
of them. The change in sign of m can be overturned due to the symmetric sum in m in 
the Fourier sums in (5.14) and (5.17). The change in sign of u, however, is a problem when 
trying to relate a term with _iSf muJ to a term with + iSf muJ due to the non-symmetric nature 
under (m,uj) — > (— m, —uj) of the integrals over w or w for all states involved in (5.14) and 
(5.17). It follows that A (i muJ T tr + Vi mu) T tr ) is not constant in r and therefore neither is 
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A (Ttr) nor A (T tr \ . Similarly, a calculation of (i mtJ Tte + VimJTte) shows that it is 

\ / rcn \ / rcn 

not zero and therefore neither \Tte) nor \Tte) are zero. 

\ / ren \ / ren 

Indeed, Graphs 5-12 numerically corroborate the above conclusions. Graphs 5 and 6 show 
that neither A (T tr ) nor A ( T tr \ are constant in r. Graphs 9-12 show that 

I* \U~-E- / A \CCH--U- 

neither \Tao) nor \Ta6) are zero. Graphs 7-8, however, seem to indicate 

\ / ren \ / ren 

that both A \ T r ^j and A (Tr<t>) might actually be constant in r. 

Note that if instead of using expressions (5.14) we use CCH's expressions for the ex- 
pectation value of the stress tensor, we encounter the same difficulty when trying to prove 
whether T t e is zero and whether AT tr is constant in r. 

In the case a = 0, since the spin-weighted spherical harmonics hXim do not depend on u>, 
we only need a change in the sign of m to relate terms with | + i^ m | 2 to terms with |-iYj m | 2 in 
Ai muJ T tr . Indeed, use of (2.22) allows us to prove that ^ m Az mii; T tr is constant in r and that 
Y^m imujTtr = ~ Ylim irwJTfr m the Schwarzschild background. A numerical investigation of 
the luminosity and components of the RSET for spin-1 in the past Unruh state is presented 
in [16]. 

The solution to this deadlock for the spin-1 case in the Kerr background consists in 
integrating over the solid angle. This allows us to relate a term with J dQ^iS?^ to a term 
with f dQ + iS? mLU , when both types of terms appear in f dQAi moJ T tr . This is in accord 
with the fact that if we integrate the conservation equation (6.4a) over the solid angle we 
immediately obtain that 

J dQAT tr = J dVLK{6) = const. (6.9) 

Indeed, we found that 

J dQA lmw T^ = - j dnA lm ^ = ^ujJ lmu , (6.10) 

where we have included the constants of normalization (3.4). We can now give simple 
expressions for the luminosity when the electromagnetic field is in the past Boulware state 
and in the past Unruh state: 

dM 



dt 
dM 



— / d^Tz mil ; (6.11a) 

71 1=1 m=l P=±l ^° 

sEEEf^ ("lb, 

1 1 « 7 D LI *J U 



dt 



1=1 m=-l P=±l 
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in agreement with Page's [37]. The former is a manifestation of the Starobinskii-Unruh 
radiation and the latter of the Hawking radiation. Since only superradiant modes are being 
included in the Starobinskii-Unruh radiation (6.11a) and the transmission coefficient Ti mu) 
is negative for these modes, there is a constant outflow of energy from the black hole when 
the field is in the past Boulware state. 

We numerically evaluated (6.11) for the case a = 0.95M. The results, compared against 
values in the literature are: 



M 2 
M 2 



dM 



dt 
dM 



dt 



B~ 



B~ 



— 4.750 * 10"' 
-5.01 * 10" 5 



(spin-1) 
(spin-0, Duffy) 



(6.12a) 
(6.12b) 



in the past Boulware state, and 
dM 



M 



M' 



dt 
dM 



dt 



u- 



u- 



-1.1714* 10' 



-1.18* 10" 



(spin-1) 
(spin-1, Page) 



(6.13a) 
(6.13b) 



in the past Unruh state. The value (6.12b) for the scalar field is calculated by Duffy [32] 
and we have calculated (6.13b) from splining Page's [37] numerical results. Both of them 
have also been calculated for a = 0.95M. 
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from that of the previous graphs in order to make more visible the region far from the horizon. 
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VII. RSET CLOSE TO THE HORIZON IN THE BOULWARE VACUUM 

Candelas and Deutsch [38] consider flat space-time in the presence of an accelerating 
barrier with acceleration a^ 1 . They then calculate the spin-1 RSET in the tetrad of an 
accelerating observer RO with local acceleration In the limit ^/as — > oo the vacuum 
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state above the accelerating mirror approximates the Fulling vacuum \F). The result is 



where the bars on the indices indicate RO tetrad. Expression (7.1) is equivalent to minus 
the stress-energy tensor for thermal radiation at a temperature of (27r^)~ 1 . As mentioned 
in Section IV, in the Schwarzschild space-time, analogously to (7.1) in flat space, the RSET 
close to the horizon when the field is in the Boulware vacuum diverges like minus the stress 
tensor of black body radiation at the black hole temperature. It is therefore reasonable to 
expect that if there existed a state in Kerr with the defining features that the Boulware 
vacuum possesses in Schwarzschild, then the RSET close to the horizon when the field 
were in this vacuum, would diverge like minus the stress tensor of black body radiation at 
the black hole temperature rotating with the horizon. The past Boulware vacuum is not 
invariant under (t, 0) reversal because of the existence of the Starobinskh-Unruh radiation. 
However, the stress tensor components tr and r0, which correspond to the Starobinskh- 
Unruh radiation, are expected (from Section VI) to have a divergence of one lower leading 
order than that of the diagonal components as the horizon is approached. It is with this 
understanding that we say that a state is isotropic at the horizon and that, in particular, 
the past Boulware vacuum might be isotropic. It is obvious that to next order in A the past 



invariant under (t, <j)) reversal either. 

CCH claim that the RSET of the electromagnetic field in the past Boulware vacuum close 
to the horizon differs from that of minus the stress-energy tensor of a thermal distribution 
rotating at the angular velocity of a Carter observer by a factor which is a function of 9. 
In the present section, we will show that CCH's result is due to a flawed assumption in the 
asymptotic behaviour of the SWSH. We will show this by re-calculating their result using 
the assumptions we believe they used. The numerical results back up the fact that the 
mentioned RSET is (minus) thermal at the horizon. 

We also mentioned in Section IV that Frolov and Thorne claim that close to the horizon 
ZAMOs measure a thermal stress tensor which is rigidly rotating with the horizon when the 

/ A y FT—B 

field is in the \FT) state. That is, they argue that (T^ u ) (where \E>) is an unspecified 




(7.1) 




might be since \CCH ) is not 



ren 
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Boulware-type state) is thermal close to the horizon, and isotropic in the frame of a RRO. 
Duffy, in turn, shows that close to the horizon RROs measure a thermal state which is 
rigidly rotating with the horizon when the field is in the \Hm ) state in the Kerr space-time 
modified with the introduction of a mirror. That is to say, close to the horizon ( T^ v ) 

\ / rcn 

is thermal and isotropic in the frame of a RRO. Finally, as mentioned above, CCH claim 

/ „ \ CCH--B- 

that ( T^y ) at the horizon differs by a factor from a thermal distribution isotropic 



in the Carter tetrad. Of course, the angular velocity at the horizon of a Carter observer, a 
RRO and a ZAMO is Q = Q + for them all, so that CCH's result does not actually distinguish 
between these observers. 

Ottewill and Winstanley [39] have proved that if a certain stress-energy tensor is thermal 
and rigidly-rotating with the horizon everywhere, then it is divergent on the speed-of-light 
surface in the Boyer-Lindquist co-ordinates, which are regular on this surface. This implies 

/ „ xCCH~-B~ 

that if ( T^y \ were thermal and rigidly- rotating with the horizon everywhere then 

\ / rcn 

the state \CCH~) would have to be irregular on the speed-of-light surface. In the present 
section we will numerically investigate the rate of rotation of the thermal distribution in 
question. 

The stress-energy tensor of a spin-1 thermal distribution at the Hawking temperature 
rigidly rotating with the horizon is given by 



1 1 TV 2 
^(th.RR)^ _ lL1 n 

v ~ 45 

where 



X P X P 



(7.2) 



T = — %= (7.3) 

is the local temperature. Note that this stress-energy tensor is obviously isotropic in the 
frame of a RRO, but it is not in the rigidly-rotating co-ordinate system {t + , r, 9, 0+}, where 
+ = — Q + t and t + = t, which is not adapted to a RRO. 

In primed co-ordinates, which are adapted to a RRO, the rigidly-rotating thermal stress 
tensor becomes 

T (th,RRV = H(r+-r_) 4 1 / 1 1 1\ 

v 28.3 2 -5tt 2 A 2 ^ 2 ^ g V '3'3'3j 1 j 

in the Kerr space-time. 

CCH calculate an expression for the RSET close to the horizon when the electromagnetic 
field is in the past Boulware state. They make the assumption that the RSET close to the 



39 



horizon when the field is in this state is more irregular than when it is in the \CCH~) state 
and therefore approximate 

A \B~ I „ \ B~ / „ v CCH- / „ \ B~ J „ v CCH- 

T» v ) ~{T» V ) -(t» v ) = (tM -It\) (r^r+) (7.5) 

/ ren \ / rcn \ / ren \ / \ / 

One can then use the expressions for the expectation value of the stress tensor, and it is clear 
that only the 'up' modes are involved in the calculation. Their result, when the components 
of the stress tensor are put in the Carter orthonormal tetrad is: 

/ren 7T 2 A 2 Zj e 2 ™' K - 1 & V 3' 3' 3, 

1 llfr, -r) 4 1 / 1 1 l x 



r+ 2 8 -3 2 -5tt 2 A 2 S(2Mr + ) diag V 3' 3' 3, 
where the hats on the indices indicate adaptation to the Carter orthonormal tetrad. This 
expression and the expected result, minus (7.4), differ in a factor of r + (2Mr + )/S. We 
proceed to reproduce CCH's expression to explain this disagreement. 

We believe that CCH followed Candelas [6] method for spin-0 to obtain asymptotic 
expansions for the radial solutions for spin-1 close to the horizon. This is the method that 
we develop for spin-1 in Appendix B. Armed with the asymptotics of that appendix, we 
can proceed to calculate the different components of the stress-energy tensor. In order to 
do that, we are first going to separately calculate the asymptotic expressions for the various 
terms that occur in the classical stress-energy tensor (2.23). 

As mentioned in Appendix B, for the asymptotic behaviour we are seeking here we 
can replace the spin-weighted spheroidal harmonics hZi mu} by the spin-weighted spherical 
harmonics hXim- Note that this implies that the result of this asymptotic analysis is the 
same whether it is equations (5.14) and (5.17) or CCH's expressions for the expectation 
value of the stress tensor that are used. We can make use of (2.21), which immediately leads 
to 

i 

^ imu^h 'imw^T ^° r -> r +) when h^ti (7.7) 

m=— I 

The asymptotic calculation of the term |z mc j0o P | 2 requires a more careful treatment. We 
observe in Appendix B that the large-/ modes dominate the Fourier series for the 'up' radial 
solution close to the horizon. We can therefore replace Yl^o with J °° dl and use the fact 
that of the two independent variables uj and m, hRfmui depends only on Co in the limit 
(I — > +oo,r — > r + ), whereas hXim depends only on m. From (2.15) and using equations 
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(B8), (Bll), (2.22), (2.21) and other properties of the spin-weighted spheroidal harmonics 
that can be found in [181 we then obtain that 



Lm,P 



2%S 



d// 3 |iV^| 2 



[r — > r A 



We substitute (B9), (BIO), and (B12) in the above equation and approximate \Bi muJ ~ 
which is valid in the limit of large /. The next integral, found in [40], is needed: 



(7.8) 
I 2 - 



d^(2^/ 2) = ^ 2 (i + 0|rN)r 



We finally obtain 



up 1 2 



3 ■ 2 A x 2 



Mr+u) Itftl 



r — > Tj 



(7.9) 



(7.10) 



6vr 2 SA 2 

l,m,F 

The other terms in the expression for the stress-energy tensor can be obtained in a similar 



manner, but they are easier to calculate. We will therefore only give the final results: 



E 

l,m,P 

E 

Lm,P 



i up 
lmu><P—l 



2 2Mr+u 91 



j^p I 



3tt 2 A 3 

Mr + tu |9l| 
6vr 2 S 2 A 



(7.11a) 
(7.11b) 



We can now use equations (7.7), (7.10) and (7.11) together with the quantum expressions 
(7.5) and the expectation value of the stress tensor to reproduce equation 3.7 in CCH. We 
obtain 



B~ 
ren 



B~ —C'CH~ 



-8M 3 7 



duouj(uj 2 + k 2 



X 



(7.12) 



3vr 2 A 2 E 2 

/-3(r 2 +a 2 ) -a 2 sin 2 6 4a sin 2 6{r 2 + + a 2 ) \ 
SO 

OS 

y -4a (r 2 + a 2 ) + 3a 2 sin 2 9 J 

in Boyer-Lindquist co-ordinates. This is exactly equation 3.7 in CCH except for the fact 
that (7.12) contains a factor r\ instead of a r + in CCH. We believe that the discrepancy 
is due to a typographical error in CCH since otherwise the stress-energy tensor would not 
have the correct units. We have also checked that equation (7.12), when the tensor indices 
are adapted to the Carter orthonormal tetrad, produces the result (7.6) above. Again, 
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the discrepancy with respect to (7.6) is only in the power of r + . It seems that, despite 
the dicrepancy in the power of r + , this is the method that CCH used to calculate their 
expression (7.6). However, as we point out in Appendix B, this asymptotic analysis is only 
valid when both uj and m are kept bounded since otherwise we would not be able to replace 
the spin-weighted spheroidal harmonics by the spin-weighted spherical harmonics. In the 
analysis we have just carried out uj and m do not both remain bounded in general. The 
only points in the Kerr space-time where both remain bounded are the points along the axis 
9 = or 7r since, there, the Newman- Penrose scalars imw^h are only non-zero for m — ±1, 
and thus m is bounded. The frequency uj is then also kept bounded because the factor in 
the integrand diminishes exponentially with uj and thus the contribution is only important 
when uj is bounded. Equations (7.12) and (7.6) are therefore only valid at the axis. An 
asymptotic behaviour of the 'up' radial solutions uniform both in I and uj is required. 

Another issue is the fact that the state \CCH~) has been used in (7.5) as a Hartle- 
Hawking state, regular on both the past and future horizons. We know from Kay and Wald's 
work that there exists no such state on the Kerr space-time satisfying its isommetries. Since 
\CCH~) is not invariant under (t,(j>) reversal it is not covered by Kay and Wald's result 
and thus it might be regular on both 7i~ and 7i + . We saw that Ottewill and Winstanley 
[31] argued that in the scalar case this state is irregular on 7i~ and regular on 7i + . Even if 
that were also the case for spin-1, using \CCH~) in the preceding calculation could still be 
acceptable if the divergence of \CCH~ ) close to r + is of a smaller order than that of \B~ ). 
In the Schwarzschild background Candelas has shown that the Unruh state is irregular on 
H~ , regular on 7i + and that the order of its divergence close to r + is smaller than that of 
the Boulware state. It is therefore reasonable to expect that the order of the divergence of 
\CCH~ ) close to r + in the Kerr background is smaller than that of the past Boulware state. 
Indeed, our numerical data indicate that the approximation in (7.5) is correct. 

Graphs 13-17 show that the RSET when the field is in the past Boulware vacuum ap- 
proaches a thermal distribution rotating with the horizon rather than CCH's result (7.12). 
The red lines in the graphs correspond to the thermal stress tensor (7.2) rotating with 
the horizon located at r = r + ~ 1.3122. The black lines are also located at r = r + 
and correspond to CCH's result (7.12). It can be seen in the graphs that as r becomes 



rather than CCH's corrected equation (7.12) (black line). At the poles, however, it can 




approaches the thermal stress tensor (7.2) (red line) 
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be straight-forwardly checked analytically that the two coincide, as expected. Only for the 
rr-component, which is the only component that diverges like 0(A -3 ) close to the horizon, 
we did not seem to be able to obtain a clear plot, which we do not include. 

Within the range of r considered in Graphs 13-17 (except 16) for the difference between 
the states \CCH~) and \B~) of the various expectation values, the corresponding plots 
for the difference between the states \U~) and \B~) are identical. This is the expected 
behaviour since for small radius r the 'up' modes dominate in these RSETs. Graph 16 
includes the two differences for the ^-component of the stress-energy tensor up to a value 
of r large enough so that the two differences become clearly distinct. 

In following with the notation used in (2.2) and the one used so far for tensor components 
in Boyer-Lindquist co-ordinates, we use the obvious notation of '(a/3)-component' to refer 
to the stress-energy tensor component T^e^^e^ in the tetrad of a stationary observer. 
Since the angular velocities of a RRO, ZAMO and Carter observer all equal Q + at the 
horizon, each one of the diagonal components of a stress tensor for a thermal distribution 
will be the same in any of the three tetrads adapted to these observers. The (r0)-component 
will also be the same in any of the three tetrads since the tetrad vectors ei r ) and e<Q\ do 
not depend on the rate of rotation. The (t0)-component, however, vanishes to leading order 
for the radial functions as r — > r + . To the next leading order for the radial functions this 
component does depend on the rate of rotation of the stationary observer that the tetrad 

/ „ \ CCH--B- 

is adapted to. Graphs 18-19 for \T t , ) show that the rate of rotation of the 

\ / rcn 

thermal distribution approaches, to next order in A, that of a RRO, rather than that of a 

ZAMO or a Carter observer. This result tallies with Duffy [32] 's results for the spin-0 case 

in the Kerr space-time modified with a mirror when the field is in the \Hj^) state. He also 

/ " \ U ~~ B ~ 

numerically shows that ( T^ u j is, close to the horizon and for the scalar field, thermal 

\ / ren 

and rotating at the rate of a RRO to 0(A) in the angular frequency. We calculated and 

/ A \U-~B- I „ xCCH~-B- 

plotted (Tt + <f> + ) and it fully coincided with \ T t+( p + ) in the region of Graphs 

\ / rcn \ / rcn 

18-19, which is why we do not include them. We conclude that the rate of rotation close to 
the horizon for the difference between the states \U~) and \B~) is also that of a RRO, in 
agreement with Duffy's results. 

An alternative technique for investigating what is the rate of rotation of the thermal 
distribution at the horizon is as follows. We find what is the frequency Q = ^zefo °f 
rotation of the tetrad frame (2.2) such that T^) = 0, where the term ZEFO stands for zero 
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energy flux observer. The answer is 

— 2(7 

fi Z EFO = ; (7-13) 

where 

B = g^Tu - gttTtf, (7.14) 
C = gt<t>Ttt — gttTty 

i „ \CCH~-B- 

We then plot I^zefo where is replaced by ( T M „ ) in (7.14). This plot is compared 

\ / ren 

against that of the angular velocities of a RRO, ZAMO and Carter observer in Figure 20. 

/* \ U '~ B ~ 

We also plotted ^zefo where is replaced by ( ) and it fully coincided with the 



ren 

CCH--B' 



corresponding one for I ) in the region of Figure 20. 

\ / ren 

Graphs 21 show the behaviour of the various modes as the horizon is approached. The 
behaviour of the 'up' modes close to the horizon are explained by the horizon asymptotics 
developed in Appendix B. 
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(cyan). 




(th, Carter) 



(yellow). The light blue and brown surfaces correspond to the speed-of-light and the static limit 
surfaces respectively. 
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FIG. 20: Plots of (f^zEFO — ^+) and (^zefo — &+) /A (dark surfaces) where 7]^ is replaced by 

A \CCH--B~ 

T„ v ) in (7.13), together with the corresponding plots with the angular velocities of a 



RRO (red), ZAMO (magenta) and Carter observer (yellow). 
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FIG. 21: The sum over Z has not been performed. For each value of I the sum Yum=—l nas been 

/ „ \CCH--U- 

performed. In the case of ( Tee / the low-/ modes clearly dominate close to the horizon. 



On the other hand, the high-Z modes dominate close to the horizon in the case of ( 



U--B- 
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APPENDIX A: RADIAL FUNCTION 



This Appendix is valid in the Kerr-Newman background with the understanding that 
the background is considered fixed and the NP Maxwell scalars include the electromagnetic 
field associated to the background as well as an electromagnetic perturbation. The quantity 
Q represents the charge of the Kerr-Newman black hole as viewed from infinity; in this 
Appendix, the quantities A and r± take the following values: A = r 2 — 2Mr + a 2 + Q 2 and 
r± — M ± y/ M 2 — a 2 — Q 2 . The Kerr background is recovered in the limit Q = 0. 

The potential in the radial Teukolsky equation is a long-range potential. Our numerical 
calculations followed Detweiler [41], who derived from the radial Teukolsky equation another 
differential equation with a real, short-range potential. We numerically solved this equation 
and obtained its solution Xi muj and its derivative. We then find from these the radial 
functions +iRi miJ , -\Rimw an d its derivative with expressions in [41]. The independent 
variable in Detweiler's differential equation is the tortoise co-ordinate r* defined by dr^/dr = 
(r 2 + a 2 )/A. Detweiler's potential tends to — uj 2 + 0(r~ 2 ) at infinity (r* — > +oo) and to 
— Co 2 at the horizon (r* — > — oo), where u = uj — mQ + . We can therefore define two sets of 
solutions with the following asymptotic behaviours: 



When the behaviour of the solution modes in terms of the time t and the angle <fi is included, 
we can find the asymptotic behaviour of the solution modes in terms of the advanced (v = 




+oo) 



+oo) 



(Alb) 



(Ala) 
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t + r*) and retarded (u = £ — r*) time co-ordinates: 



x^in +imcj>—iujt 



Q—iu)v+im<f> 



at X 



^in ^—iLdu+im<j> ^ 



at 7Y" 



5 Le" ! ™ + ""^ + at 



(A2a) 



\^up +im(h—iu}t 



5 



up 

Imuj 



at X 



-vuju-\-vm4> q± 2"+ 



(A2b) 



A?rL> e 



at 

at W+ 



Equation (A2a) represents a wave emerging from X~, being partially scattered back to X + 
and partially transmitted through to 7i + . Similarly, (A2b) represents a wave emerging from 
T~C~, being partially scattered back to 7i + and partially transmitted through to X + . 



Both sets of modes are eigenf unctions of the hamiltonians Hk y = iky and 



ik 



3d? 



with eigenvalues to and uj respectively: 



tt v» +im4>—iiot 



V"» +im<f>—iu)t 



£t v „+imd>—ibjt ~. v „+im<h—iujt 



(A3) 



where the symbol • indicates either 'in' or 'up'. We will restrict the definition of the 'in' 
and 'up' modes to those modes with positive uj and positive uj respectively. It then follows 
that the 'in' and 'up' modes are positive frequency with respect to the Killing vectors k j 
and kjtf respectively. 

The asymptotic behaviour of solutions hRp^ V7i = 0, ±1/2, ±1, ±3/2, ±2 of the radial 
Teukolsky equation is: 



in,tra a —h„—iwr* 



hR 



pm,tra a 



irnw 



m,mc i -iur* i r>in,rei —l-2h „+iu>r* 



Imuj 



hRlm 



h^lmu e + hR lmuJ A e 



tmu) 



hR,, 



up,tra -l-2h. 



Imuj 



r 



+oo) 

r+) 
±oo) 



(A4a) 



(A4b) 
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From (Al) and expressions in [41] we can find the asymptotic coefficients of _i-R* mu; from 



those of X' muJ : 



j->in,rcf 
-l U lmu> 


4a; 2 


7->in,tra 


—sgn(u) \oj\i 


pin.inc *\ n 
T-)Up,rcf 


l-Blmu) 


Din.inc r^i n 
pup,tra 


(r 2 + a 2 ) 1 ^* 
M(r 2 +a 2 ) 1 / 2 


oup.inc /.up 

pin,inc 


i 


pup,inc niip 
pup,inc 


l^+l 

-2 X /2( r 2 + Q 2)l/2£, 


2 3 / 2 |a;| 





(A5) 



In the calculation of _ii?J°^ a /_i-Rj^ c in (A5) we needed an extra term in the asymptotic 
expansion of the ingoing part (the outgoing part is simply obtained by complex conjuga- 
tion since Detweiler's potential is real) of Xi mu} close to the horizon. By introducing the 
asymptotic expansion 

-jffi" = [1 + ai (r - r + ) + 0((r - r + ) 2 )] e"^* (A6) 

Imuj 

in Detweiler's differential equation and performing a Taylor series expansion around r + of 
Detweiler's potential, we find from the second order term that 



0L\ = 

AMa 2 r + -Q 2 (a 2 - r 2 ) a 2 + Q 2 4amr + Co 2(r + - M) 2 V 



2m* 



-lAzmu; 



(A7) 



{r\ + a 2 ) 2 r\ + v 2 r + — r_ (r+ + v 2 ) 2 

where v 2 = a 2 — am/uj and r\ = (\Bi muJ — _i\i muJ )/ (2u 2 ). In the calculation of _i-R|^™ c , an 
extra term is also needed in the asymptotic expansion of the ingoing part of Xi muJ for large 
r: 



X-lmui 

with 



l + ^ + 0(r" 2 ) 
r 



e~ iuir * (A* 



a i-ikmw + 2aurn)i 

Pl ~ 2^ (A9) 

After obtaining the asymptotic coefficients of _iRJ mw from those of X' mtJ , we just need 

to derive those of + \R' muJ to complete the asymptotic picture of the solutions to the radial 

Teukolsky equation for spin-1. This is achieved by using the asymptotic behaviour in (A4) 

together with a transformation in [41] that relates + \R to _i-R; the result is: 



j->in,inc j->in,ref T-)Up,tra rj2 

+^ n imu _ _0 3 ij 2 - + lU lmuj _ +^ n imu _ l-^lmu 

^in.inc ' R m,rc{ fiW,** 2uJ 2 

1 -l u Imu _1 !mu 

7-)in,tra j-)Up,ref T->up,inc ■ r>2 

+l n irnoj _ +l U lmuj _ _2 3 J{ fftV + lJX lmm _ ~ l ^ a lmui 

D in,tra D up,ref D up,inc 9 ts m 



(A10) 
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It is clear from their asymptotic behaviour that neither h-RJ^ nor hRfmui satisfy the 
symmetry (2.10a). As a matter of fact, we obtained that, under this symmetry, the functions 
hR'mu transform to the radial funcions that are derived from the solution X*^ and its 
derivative. Indeed, we construct first a new radial function _iR' lrnu] derived from X'^ in the 
same manner that -iR° mui is derived from X' muJ , using an expression in [41]. It can then be 
checked using equations (A5) and (A10) that the relations 

i .* (AU) 

are satisfied, where +iR' muJ is calculated by applying to _ 1 ^J' mw the transformation in [41] 
that relates + ±R to and its derivative. Renaming ±iRj^ luJ and ii-R^L by ±iR?^ and 

pdown 



-i-Rfmw 1 respectively, we have the following two sets of modes: 



±1 iCl = (2 lJ B^) ±1 A^ 1 TlJ R)- ~ (A12a) 



(Al2b) 



MjCr-^ + iiC^-V- (r - +oo) 
±1 RfZ n = (^B^A^^RfZ ~ (A12c) 
' ± 1 ifiT toc A^e-* 5 " + ± i<7' rcf e + ^* (r - r + ) 

±ijR down,t r a r _ le _^ (r ^ +Qo) 



(A12d) 



where 



jr-,out,inc/rcf/tra /o D T->in,inc/ref/tra* 7~)down,inc/ref/tra /o D j-)Up,inc/rcf/tra* 

(A13) 

Note that the factor (2iS; mw ) ±1 is needed so that the Teukolsky-Starobinskh identities are 
satisfied. Similarly, since the radial modes ±i-R°^ and ±i-R^ n are obtained from A^ and 
A^* respectively, we rename the latter as 

= ^ (A14a) 

+iujr * + A^e~ iujr * (r -> +oo) 

= ~ ^ lmuj (Ai4b) 
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with 



a out _ a in* 
Irnui imuj' 

rjout _ nin* 

lmu> irnti" 



A 



B 



down 
Imu) 

down 
Imw 



up* 



DIP* 

Imu} 



(A15) 



We follow the same positive-frequency convention for the 'out' and 'down' modes as that 
for the 'in' and 'up' modes respectively, namely, their definition is restricted to modes with 
positive uj and positive u respectively. The asymptotic behaviour in terms of the advanced 
and retarded time co-ordinates of these two new sets of functions is 



X 



out +im(j>— itot 



4 out 



ydown +im(f>— iuit 



iujv+im(f> 

^—iuju+im<f> 

rjout „— iCou+imcf)^ 
- D lmuj c 





rjdown _— iojv+im(f> 
D lmuj e 



at X 
at X+ 
at Ti~ 
at H + 

at X~ 
at X+ 



y^downg— iuiu+im<j)+ g^- 



(Al6a) 



(A16b) 



-iCbv+im<f>j r 



at 7i~ 



Modes (Al6a) describe a wave going out to X + whereas modes (Al6b) describe a wave going 
down HA . 

It is easy to check that the 'out' ['down'] NP scalars are related to the 'in' ['up'] NP scalars 
under the symmetry transformation (t, 0) — > (—t, — 0) in the manner: 





6 in/up 








/in/up 








/in/up 


Intuit 


Al 



.l)m+l A -l p -2 



, out / down 
I— m— uj(P-\-1 



■(-!)' 



m+l ,out/down 
1—m—uj'rQ 



> under (t, 



. 1)m+ i ApVm _ 



i out /down 



(A17) 



where we have used the radial symmetry (2.10b) and the angular symmetry (2.11c). 

Finally, it can be checked that the radial Teukolsky equation admits the following wron- 
skian: 



(A18) 



where C" are real constants. It is useful to note the following two relations satisfied by the 
wronskians of the 'in' and 'up' solutions once the normalization constants (3.4) have been 
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included: 

W[+iiC. -liffl = -i*t*J (A19a) 

W / [+l^"m^) = +W / [+l-R*-m-a; ) -l^Z-m-J (Al9b) 

As a short-hand notation we will denote W[ + iR' miiJ , -i-R'^J by W[ + iR, -iR*}' muJ . We give 
here two important wronskian relations, obtainable in terms of the coefficients of the solution 
to Detweiler equation by using relations (A10) and (A5): 

u^-\XLS)=u\B? m S (A20a) 
^(l-KL| 2 )^IC| 2 (A20b) 

I 1 2 

Wronskian relation (A20a) shows that for the modes such that ujuj < 0, it must be LAjJJ^J > 
1 , and thus the wave mode is reflected back with a gain of energy. This phenomenon is 
known as superradiance. Since the 'in' modes are only defined for positive u, superradiance 
occurs for these modes for negative u only. Similarly, from wronskian relation (A20b), when 
Cjuj < it must be |AJ^J 2 > 1. Therefore, the 'up' modes, which are defined for positive 
G), that experience superradiance are those for which u < 0. The condition Colo < for 
superradiance, which is the same for scalar and gravitational perturbations, clearly shows 
that this phenomenon is only possible if a ^ and therefore it only occurs if the black hole 
possesses an ergosphere. 
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APPENDIX B: ASYMPTOTICS CLOSE TO THE HORIZON 

This Appendix is valid in the Kerr-Newman background with the same re-definitions of 
quantities as in Appendix A. 

In this Appendix we are interested in finding the asymptotic behaviour close to the 
horizon of the solution to the radial Teukolsky equation. The main application of this study 
is in the calculation in Section VII of the behaviour close to the horizon of the RSET when 
the field is in the past Boulware state. We will therefore only consider the 'up' modes and, 
because of the presence of a factor that decreases exponentially with u, we will not consider 
the case of large Co. 

This study is based on one performed by Candelas [6]. Even though Candelas started 
the calculation for general spin, he soon confined it to the scalar case. It is our intention to 
complete his asymptotic calculation for general spin and only specialize to the spin-1 case 
at the end. 

Candelas performs a Taylor series expansion around r = r + of the coefficients of hRimu 
and its derivatives appearing in the radial Teukolsky equation. By keeping only the first 
order terms in the expansion and also terms that involve parameters which might become 
very large, the radial Teukolsky equation becomes: 

hRlmui = (Bl) 

where q = 2K + /[r + — r_). The parameters in (Bl) that might become very large indepen- 
dently of the limit r — > r + are h^imuj, & and Co. We keep Co bounded, which means that either 
both m and uo are bounded or else that to — ► oo and m ~ oo/Q + . We are going to restrict 
ourselves to the first possibility (i.e., m and io bounded) since it is only for this case that we 
are able to find the behaviour of the angular solutions, needed in the calculation in Section 
VII. We thus have that the only term in (Bl) that might become very large independently 
of r — > r + is the one with hXimw Since we are keeping m and u bounded, hA; maJ can only 
become large if we let I — > +oo. 

When letting / — > +oo and keeping lo and m bounded in the Teukolsky angular equation, 
all the terms in the coefficient of the angular function h£>imu(9) can be ignored except for 
h\mui an d those with a 1/ sin 9 in them. This is equivalent to setting au = in the angular 
equation. This means that in the limit I — > +oo (with m and to bounded) the angular 



/ \ d hRlmu) .11 . -i \ ^-hRlmuj 

[r-r + ) — +(h+l)- 



dr s 



dr 



hMmw - 4:iuhr _ q(q - 2ih) 
4(r — r+) 



r _ r _ 
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solution reduces to the spin-weighted spheroidal harmonics: hZimu> hXim an d that hMmu — > 



After replacing h^imu by I 2 and neglecting Aiuhr in (Bl), Candelas goes on to prove that 
the asymptotic behaviour close to the horizon of the 'up' radial functions with / — > +00 is 
given by 



hRiL -> h aix- h/2 K h+iq (2lx 1/2 ) + h b lX - h/2 I^ h+iq) (2lx 1/2 ) (I +00, r -> r+) (B2) 



which is uniformly valid in Z. The factors ^a; and ^ are the coefficients of the two indepen- 
dent solutions. It is at this point that Candelas' analysis specializes to the scalar case. We 
pursue it here for general spin. 

The asymptotic behaviour of the modified Bessel functions I v {z) and K u (z) are well 
known ( [42]). If we fix r close to the horizon in (B2) and let I — > +00, the value of the 
radial potential at that fixed value of r will go to infinity and thus it must be h-ft)^ —* 0. 
From the asymptotic behaviour for large argument of the modified Bessel functions, it is 
only possible that hR^muj in the limit I — > +00 with r fixed in (B2) if the coefficient 
hbi decreases exponentially with /. It follows from this result together with the asymptotic 
behaviour for small argument of the modified Bessel functions that in the limits r —>■ r + 
and / — > +00, while keeping Ix 1 ^ 2 finite, the second term in (B2) can be neglected with 
respect to the first one. In this last statement we have made the implicit assumption that 
if the coefficient h^i decreases for large I, then it does slower than the coefficient ^bi- This 
assumption is proved to be correct in what follows. We have from the above discussion that 



We can determine the coefficient h^i by comparison with the WKB approximation (A4b). 



(I — h) (I + h + 1) — ► I 2 . The latter expression implies that \B, 



Imw 



I 2 in the same limit. 



hRiL - h aix- h / 2 K h+iq {21X 1 ' 2 ) (I - +00, r -+ r+, lx 1 ' 2 finite) 



(B3) 




*2sin[(/i + ig)7r]r(l - h - iq) 




(B4) 



2sin[(/i + zg)7r]r(l + h + iq) 



(I — > 00, r — > r + , lx 1 / 2 — > 0) 



where 



I Q = e~^+ (4M K+ ) 



■+ 



(-4M/c_) 2 



(B5) 



Comparing this asymptotic expression with the WKB approximation (A4b) it follows that 

2 sin[(fe + iq)ir]T{l + h + iq)I? h _ % inc 
h^i ~ z 1 hK-imu 

(B6) 

pup.rcf _ r(l + h + iq) , \2h T *2l-2(h+iq) pup.inc 

- ~ y{\ -h-iq) ^ + ~ ' " hK imco 

We now specialize to the spin-1 case. Combining equations (B3) and (B6), and using the 
same normalization as the one used in the numerical results (i.e., setting -\B^™ C = 1 and 

pup,in 
+l n lmu> 



^Z nC = -^V(2^ + )),wehave 



(B7) 



^r u p x~ 1/2 K^- (2lx l/2 ) 

+lKl ™ {r+-r-)K + T{-iq) X K ^ Mx > 

(I — > +oo, r — > r + , Ix 1 / 2 finite) 

(l — >• +oo, r — > r + , Zx 1 / 2 finite) 

It is also useful to give the expressions that the 'up' radial functions (B7) adopt in this 
limit whenever the constants of normalization (3.4) are included. These expressions are, in 
compact form: 

\N-\\ h RZj ~> A h Nx- 1/2 K h+iq (2lx 1/2 ) (I -> +oo,r -> r+,/x 1/2 finite) (B8) 
for /i = ±1, where 

f -4 ,/i= +1 1 
2z I [/(r + - r_)]- fc (B9) 

and 

T*j-iq 

N = ; g (BIO) 

N /2%#+r(-z< 7 ) 

We therefore have finally found the asymptotic behaviour of the 'up' radial modes with 
I — > +oo close to the horizon with a), and both m and u, bounded. We believe that this is 
the method behind the approximation given by CCH in their Tablell. Their result, however, 
does not exactly coincide with either (B7) or (B8) (as a matter of fact, in their table there 
is a quantity p that they have not defined and it cannot be the spin coefficient as it cannot 
have a 0-dependency). 

Note that there is no reason why the 'up' radial modes (B7) or (B8) should diverge in the 
stated limits. In fact, they clearly do not in the case of helicity —1. It is only when other 
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factors (as in (2.14)) are included that the resulting expressions we deal with in Section VII 
diverge in this limit. 

We are also interested in finding the result of applying the operator £>J A on the asymp- 
totic solution (B8). It immediately follows from (B8) and (B9) that 



m\Vl(A h R\ 



up 

mix)) 



A h N{r + — rJ)x 



-h/2 



(-~ + 1 + if) K h+iq {2lx^) + lx^K' h+iq {2lx^) 
(I — »■ +oo, r — > r+, Ix 1 / 2 finite) 



(Bll) 



When using a recurrence relation ( [42]) for the modified Bessel function, expression (Bll) 
for h — +1 reduces to 

|JV3|Z>J (A+xiC.) - -AA ^+J^ i^ (J - +oo, r r + , /x 1 / 2 finite) (B12) 

In graphs 22-23 we compare the numerical solution with the analytic asymptotic approx- 
imation (B7) we have found. These graphs show that this approximation indeed tends to 
the non-approximated (numerical) solution and that as r — > r + the approximation is better 
for the higher values of I, as predicted. 



60 



0.8- 

. : 



0.6- *. 



0.4- * 




r 



FIG. 22: , I = 4 . . . 7, m = 0, u = 0.01. Correspondence between colours and modes is 

the same as in Figure 23 (dots are the numerical solution and straight lines are the approximation 
from (B7)). 



1- 














1 kmO / 






0.8- 




41 hi 

l I3t> 


10 




0.6- 






I2m0 




0.4- 




1/ A / A 


11 mO 




0.2- 


n 











1.5 


2 


2.5 3 3.5 
r 


4 



FIG. 23: Relative error U ^— ' — ' ^ — , i = 1 . . . 7, m = 0, w = 0.01. 
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